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Abstract 

Effects of pH on adsorption and removal efficiency of lonizable organic compounds (lOCs) by environmental adsorbents are 
an area of debate, because of its dual mediation towards adsorbents and adsorbate. Here, we probe the pH-dependent 
adsorption of ionizable antibiotic oxytetracycline (comprising 0TCH2^, OTCH^, OTC~, and OTC^~) onto cyclodextrin 
polymers (CDPs) with the nature of molecular recognition and pH inertness. OTCH^ commonly has high adsorption affinity, 
OTC~ exhibits moderate affinity, and the other two species have negligible affinity. These species are evidenced to 
selectively interact with structural units (e.g., CD cavity, pore channel, and network) of the polymers and thus immobilized 
onto the adsorbents to different extents. The differences in adsorption affinity and mechanisms of the species account for 
the pH-dependent adsorption of OTC. The mathematical equations are derived from the multiple linear regression (IVILR) 
analysis of quantitatively relating adsorption affinity of OTC at varying pH to adsorbent properties. A combination of the 
MLR analysis for OTC and molecular recognition of adsorption of the species illustrates the nature of the pH-dependent 
adsorption of OTC. Based on this finding, y-HP-CDP is chosen to adsorb and remove OTC at pH 5.0 and 7.0, showing high 
removal efficiency and strong resistance to the interference of coexisting components. 



Citation: Zhang Y, Cal X, Xiong W, Jiang H, Zhao H, et al. (2014} Molecular Insights into the pH-Dependent Adsorption and Removal of lonizable Antibiotic 
Oxytetracycline by Adsorbent Cyclodextrin Polymers. PLoS ONE 9(1): e86228. doi:10.1371/journal.pone.0086228 

Editor: Arlstides Docoslls, Queen's University at Kingston, Canada 

Received June 16, 2013; Accepted December 7, 2013; Published January 21, 2014 

Copyright: © 2014 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction In any medium, provided the original author and source are credited. 

Funding: This study was supported by the National Basic Research Program of China (No. 201 3CB430403, http://www.most.gov.cn), the National Natural Science 
Foundation of China (Nos. 41 171 382, 21077020, and 21137001, http://www.nsfc.gov.cn), and the Program for New Century Excellent Talents In University 
(No. NCET-1 1-0047, http://www.dost.moe.edu.cn). We thank the Key Laboratory for Ecological Environmental in Coastal Areas, State Oceanic Administration 
(No. 201105, www.nmemc.gov.cn) for the support to X Cal. The funders had no role In study design, data collection and analysis, decision to publish, or 
preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 
* E-mail: xlyuncai@dlut.edu.cn 

a Current address: Laboratory of Instrumentation and Analytical Chemistry, Dalian Institute of Chemical Physics, CAS, Dalian, China 



introduction 

lonizable organic compounds (lOCs) occupy a large fraction of 
the pre-registered REACH compounds [1]. They have one or 
more pK^ values and are present in the form of a mixture of 
ionized and unionized species. The coexistence of these species is 
pH-dependent, which complicates the removal of lOCs in aquatic 
systems [1,2]. Particularly, veterinary antibiotics (e.g., oxytetracy- 
cline (Figure 1)), one emerging lOCs, are ubiquitously detected [3] 
and have aroused serious concerns about the spread of antibiotic 
resistance genes [4,5]. 

Adsorption is a widely used approach in the industry to remove 
organic pollutants from water [6,7]. This technique is effective for 
antibiotic removal [7]. Commonly, adsorption efficiency of 
antibiotics is controlled by solution pH. It is probably due to dual 
mediation of pH toward antibiotics and adsorbents. Specifically, 
the solution pH mediates the fraction of the species by altering 
ionization degree of antibiotics. The species of many antibiotics 
(e.g., tetracyclines [8,9], sulfathiazole [10] and sulfamethazine 
[1 1]) have been proven to dilfer in adsorption affinity. Meanwhile, 
conventional natural and synthetic adsorbents themselves (e.g., 



clays [8,12], humic substances [13,14], activated carbons [15], and 
CNTs [16]) probably suffer changes in structure and/ or chemistry 
from the variation of pH. In turn, these changes may render 
alteration of adsorption sites of adsorbents. The dual mediation of 
pH complicates adsorption behavior of antibiotics, hereby making 
an impediment to the understanding of their pH-dependent 
adsorption. 

Cyclodextrin polymers (CDPs) are growingly concerned envi- 
ronmental adsorbents, because they are endowed with a unique 
property of molecular recognition [17,18]. They are widely 
applied in adsorption and removal of various organic pollutants 
[19-26]. The polymers are composed of functional unit (i.e., CD 
monomer) and cross-linking agents (e.g., epichlorohydrin (EPI)). 
CD monomer with an axial open cavity of hydrophobic character 
(Figure 1) is capable of including organics (or a moiety) in terms of 
geometric compatibility [27-29]. This inclusion is driven by 
various intermolecular interactions (e.g., hydrogen bonding, 
hydrophobic interactions, van der Waals forces, and electrostatic 
interactions) [30-32], indicating a unique property of molecular 
recognition. For example, the species of some lOCs with one or 
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two pR'a values are included by CD in difTerent patterns [33,34] . 
This uniqueness is proved to remain in the polymers and 
considered as the primary adsorption mechanism [17-19]. 

Furthermore, the cross-linking agent EPI containing two 
reactive functional groups can form bonds with CD monomer 
and/ or itself The as-prepared polymer is a mixture of different 
CD units joined by repeating glyceryl linkers, representing a three- 
dimensional polymer network [18,19,21]. As the linkages of CD 
moieties are relatively polar, the polymer is highly hydrophHic and 
readily swells in water. This character makes adsorption sites of the 
polymer (e.g., CD cavity, pore channel, and network) easily 
accessible [18,19,21]. Additionally, the polymers are stable in both 
acidic and alkaline solutions [21], showing poor susceptibility to 
pH variations. It is due to less susceptibility of CD (pAj, = ca. 1 2) 
and isopropyl alcohol transformed by EPI (pX, = ca. 17.1) to 
ionization under general conditions. 

With these advantages, CDPs are suitable as model adsorbent 
for illustrating adsorption behavior of lOCs at varying pH. 
Therefore, this study was aimed to probe the pH-dependent 
adsorption of oxytetracycline (OTC, as model compound) onto 
CDPs and thus to select CDPs for removing the antibiotic. To this 
end, we investigated inclusion complexation (i.e., binding 
constants and complexation ratios) of CD and OTC through a 
spectroscopic titration technique and address CD inclusion 
mechanisms for the species on the basis of structiu-al character- 
ization and molecular simulation of the complexes. We conducted 
batch adsorption experiments of OTC onto CDPs for measuring 
adsorption aflfmity of OTC over pH. Mathematical equations for 
relating adsorption affinity of OTC at varying pH (or the species) 
to adsorbent properties were derived to reveal the nature of pH- 
dependent adsorption. The equations were used to select CDPs to 
remove OTC at different pH from water. 

Materials and Methods 

Materials 

Oxytetracycline dihydrate (OTC, purity >95%) was obtained 
from AF Pharma LLC (US). Information on four cyclodextrin 
monomers (i.e., P-CD, RMCD, HPCD and y-CD) were listed in 
Table SI. fi-Cyclodextrin was recrystallized twice in deionized 
water and dried under vacuum. Other cyclodextrins were used as 
received. Other chemicals used were of analytical reagent grade. 

Adsorbent CDPs were synthesized in our previous work [19] 
where the ratio of cyclodextrin to cross-linking agent EPI was set 
at 1:5. They comprised four onefold polymers (i.e., P-CDP, y- 
CDP, HP-CDP and RM-CDP) with single CD as complex and 
three composite polymers (i.e., P-y-CDP, fi-HP-CDP and y-HP- 
CDP) with an equimolar mixture of two CDs as double complexes. 
Their physiochemical properties (including CD content, cross- 
linking degree, swelling ratio, particle size, pore size, surface area, 
and pore volume) were characterized and compiled in Table 1 in 
reference 19 [19]. The parameter CD content of the as-prepared 
polymers varied from 38.62% to 71.91%, depending on type and 
number of the complex used. The corresponding cross-linking 
degree, representing the mean amounts of EPI per CD unit, was in 
the range of 19.82 to 4.57. 

Inclusion complexation between CD and OTC at varying 
pH 

The solutions of OTC and CD were adjusted to five pH values 
ranging from 4.17 to 9.30, by adding appropriate amounts of 1- 
mol/L HCl or NaOH. The concentration of OTC was 20 |imol/ 
L. A series of concentrations of CD were set, e.g., 0 to 2000 |imol/ 
L for P-CD, 0 to 5200 (imol/L for RMCD and HPCD, and 0 to 



Table 1. 1:1 inclusion constants (Kq, L/mol) of CD with the 
species^. 





P-CD 


RMCD 


HPCD 


y-CD 


0TCH2* 


712 


2768 


_b 




OTCH* 


554 


54 


259 


128 


OTC" 




408 


891 




OTC^" 


200 






1098 



^The Kc of the species was obtained on the basis of the spectroscopic 
measurement of apparent Kc values of OTC and LINGO optimization; 
"^approximately zero. 
doi:l 0.1 371 /journal.pone.0086228.t001 



6000 |J.mol/L for y-CD. The solutions were incubated at room 
temperature (ca. 15°C) and protected from light. Absorption 
spectra of the solutions were recorded by a UV-2300 spectropho- 
tometer (Techcomp Limited, Shanghai, China). Peak wavelengths 
varied from 353 nm to 372 nm, depending on pH. Data were 
used to obtain inclusion constant {Kq, L/mol) and inclusion ratio n 
(namely a ratio of CD to OTC) in terms of the extended Benesi- 
HUdebrand equation (Equation 1) [35]. The parameter n is 
commonly 1/2, 1 or 2, which should yield the highest values of 
fitting degree of the equation. As n is determined, the parameter 
R'c can be calculated as a ratio of the intercept to the slope of the 
corresponding equation (Equation 1). 



Ao-A [/jgAgL [I],Af.LKc[MI" ^ ' 

where Aq and A are the absorbance of OTC and the mixture at 
peak wavelength, respectively; [1]q and [A4\q are the initial 
concentrations of OTC and CD, respectively; As is the difference 
of the molar absorptivities of complexed and free OTC; L is 
optical distance; and n is inclusion ratio. 

Moreover, inclusion complexes of CD with OTC were prepared 
by the grinding method [36]. Briefly, equivalent molar mixes of 
OTC and CD were grinded for 1 h with addition of three-fold 
deionized water. The complexes were lyophUized and character- 
ized by FTIR (Shimadzu IRprestige-21, Japan) and NMR (Varian 
INOVA 400 MHz, USA). The complexes were prepared as KBr 
disks for FTIR measurement. All NMR spectra of the complexes 
were recorded in D2O with the chemical shift of HOD as 
reference. 

Adsorption of OTC onto cyclodextrin polymers at varying 
pH 

Batch adsorption experiments were conducted in 25-mL glass 
centrifuge tubes and at 25°C, with a soUd-liquid ratio of 1:50. 
Concentrations of OTC varied from 1 to 30 mg/L. The solutions 
were adjusted to five pH values (varying from 4.65 to 10.11) by 
adding 1-mol/L HCl or NaOH. Appropriate amounts of 3-mol/L 
NaCl were added to all the solutions to obtain 0.02 mol/L Na"*" or 
CI , in order to minimize interference from changes of ionic 
strength upon the pH adjustment. Then the solutions were kept in 
the dark and shaken by a rotary shaker (180 rpm) for 30 min. The 
30-min contact time was sufficient to ensure adsorption equilib- 
rium on the basis of the adsorption kinetic experiments (Text SI). 
The samples were centrifuged at 4000 rpm for 15 min. A drop of 
6-mol/L HCl was added to the supernatant (acidified to pH ca. 
2.5), in order to prevent OTC decomposition (Text S2). The 
supernatant was filtered through a 0.45-|im MUlipore membrane. 
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OTC in the filter was measured on an L-2000 HPLC (Hitachi, 
Japan). Tlie stationary phase was a HypersU ODS C-18 column 
(4.6x250 mm, 5 |im). The mobile phase comprised 22% aceto- 
nitrUe and 78% NaHjPOi-HgPO* bufier solution (pH 2.30). The 
flow rate was set at 1.0 mL/min. The column oven was 
maintained at 25°C. The injection volume was 10 [iL. The 
detection wavelength was set at 355 nm. OTC was eluted at about 
5 min. Data was processed with Langmuir and FreuncUich models 
to obtain adsorption constants (Text SI). And the adsorption 
distribution coefiicients [K^, L/Kg) were also calculated (Text SI). 

Additional adsorption experiments were carried out at three pH 
values (i.e., 5.0, 7.0 and 9.0) for all the polymers, to confirm 
loading of OTC at varying pH onto the adsorbents. The initial 
concentration of OTC was set as 20 mg/L. Once approaching 
adsorption equilibrium, the polymers were centrifuged and rinsed 
with 1 0-mL deionized water. Then, the polymers were lyophUized 
for FTIR characterization. 

Computational simulation of inclusion complexation of 
CD and OTC 

The atomic coordinates of P-CD (refcode POBRON) and y-CD 
(refcode CIWMIEIO) were selected from Cambridge Structural 
Database (CSD). Starting geometries of RMCD and HPCD 
molecules were built on the basis of the structure of P-CD using 
Chem3D Ultra software (version 8.0.3). RMCD was an optional 
derivative of P-CD with methylation of all OH groups at C2, C3 
and C6 positions (Figure 1). HPCD was the hydropropylated 
derivative at C2 position (Figure 1) with substitution degree of 6.4 
(Table SI). Starting geometries of the species of OTC were also 
built using the same software. Energy minimization of all the 
molecules was separately performed using PM3 semi-empirical 
method of Gaussian 09 software [37]. Based on the optimized 
geometries, molecular docking between CD and the species of 
OTC was carried out on CDOCKER module of Discovery Studio 
(DS, version 2.5) with the CHARMm force field. Ten most stable 
poses are obtained by the molecular docking technology. Among 
them, the pose of CD-OTC in which inclusion sites and binding 
orientation are most consistent with those of the complex 
evidenced by FTIR and NMR characterization would be chosen 
for further optimization. The selected complex structures were put 
into Solvation module of DS with the CHARMM force field to 
simulate the real water environment. The obtained structure then 
was optimized using Minimization and Minimization (QM-MM) 
protocols of DS. For the resulting structures, intermolecular 
interactions involved in the inclusion were discriminated using the 
Multiwfn 2.4 package developed by Lu and Chen [38]. In 
addition, energies of van der Waals forces and electrostatic 
interactions were recorded from the docking results of DS with the 
consideration of solvent contributions. 

Computational simulation of interactions of EPI 
oligomers with OTC 

The polymer actually is a mixture of the materials comprising 
various ratios of CD to EPI, since the polymerization is rapid and 
relatively directionless. Due to the complexity of the cross-linking 
reaction, the structure of cyclodextrin polymers remains ambig- 
uous. Here, we constructed some ohgomers of EPI (e.g., ring- 
opened, tricyclo-, pentacyclo- and heptacyclo-ones) that are 
essential units of network and pore channel of the polymers. 
These oligomers may be indicative of structural units of the 
polymers. And hence, they were used as structure mode to bind 
the species of OTC. Their geometries were optimized with the 



aforesaid method for revealing molecular interactions between the 
polymer and guest compound in a qualitative way. 

Data analysis 

All experiments were performed in triplicate and the resulting 
means were used in the following model analysis (Equation 2). 
Notably, the stability of OTC was ensured under all experimental 
circumstances (Text S2). The antibiotic comprises four species (i.e., 
OTCHg^ OTCH-, OTC" and OTC^") with three pi; (i.e., 
3.27, 7.32 and 9.11). Apparent interaction constants of OTC (i.e., 
CD inclusion and adsorption affmity) are considered as overall 
contribution of the four species (Equation 2). In Equation 2, the 
fractions of the species can be calculated by substituting pH and 
pA'a into Equations S1-S4 (Text S3). To reduce system error and 
randomness in the pH designation, we used LINGO 9.0 
optimization software (Lingo System Inc., USA) to obtain the 
interaction constants of the species. The model LINGO 1 was 
chosen and all coefiicients were constrained to be &0 (Text S3). 
The interaction constants were non-weighted. The application of 
the LINGO software would be tested by comparing the apparent 
interaction constants calculated on the basis of Equation 2 and 
obtained from the respective experiments. 

K = K+a+ +K^!i^ +K-a- +K^-a^- (2) 

where A" is the apparent interaction constant of OTC; A^, , IC~ 
and are the interaction constants of OTCH2 , OTCH , 

OTC~ and OTC^~, respectively; and a""", ot*, oT and are 
corresponding mass fractions of the species. 

Results and Discussion 

Inclusion of CD with the species of OTC 

Addition of CD caused obvious changes in the UV-vis spectrum 
of OTC (Figure SI), evidencing the formation of inclusion 
complex. Comparison of regression coefficients [R') derived from 
Equation 1 showed the complex had an inclusion ratio of 1:1 
(Table S2). The corresponding inclusion constant (i.e., iTc) was 
compiled in Figure 2 and Table S2, relating to both CD type and 
solution pH. For instance, the inclusion capacity of P-CD with 
OTC, initially increased with pH increasing, approached the 




~i ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' r~ 

3456789 10 



Figure 2. pH-dependent inclusion constants of OTC. The scatter 
represents the experimental values at certain pH and the line represents 
the calculated ones at varying pH on the basis of Equation 2. 
doi:10.1371/journal.pone.0086228.g002 
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maximum at pH 4.75, and subsequently reduced. An inverse 
profile of inclusion capacity over pH, however, was observed for 
RMCD with the minimum at pH 6.36. 

The Kq value of the species (Table 1) was obtained through the 
LINGO optimization of apparent Kq values for each CD at 
varying pH. A good linear relationship (/?^S0.74 and P<0.05) was 
observed between the apparent Kc derived from spectroscopic 
titration experiments and calculated on the basis of Equation 2 
into which the Kq and fraction of the species were substituted 
(Table S3). This relationship evidenced that the apparent Kq of 
OTC at varying pH (lines in Figure 2) could be obtained by 
Equation 2. Comparison of the values (Table 1) demonstrated 
that molecular recognition of the species was related to cavity size 
[39] (e.g., P-CD versus y-CD) and structural modification (P-CD 
versus HPCD/RMCD) of CD (Figure 1). It is due to the fact that 
the species differing in geometry (Figure 1) theoretically match the 
cavity of CD in diflFerent patterns. Overall, the differences in the 
Kq values of the species account for the pH-dependent inclusion of 
OTC, taking into account the nature of the pH inertness of CD. 

Inclusion mechanisms of CD with the species of OTC 

Both FTIR and NMR characterization also convinced inclusion 
complexation of CD (e.g., P-CD) with the antibiotic. In the FTIR 
spectrum of the complex (Figure S2), stretching vibration peaks of 
O-H (3488 cm ') and dimethylamino groups (3372 cm ') of 
OTC turned to a strong and sharp peak at 3370 cm~'. In the 
NMR spectrum of the complex (Table S4), chemical shift of the 
protons at C4 and C6 of the OTC molecule transported from 
<52.686 and <51.601 to 52.698 and 51.635, respectively. And the 
concomitants were both downshield shift of H 1 -H6 protons of P- 
CD and reduction of splitting numbers of H3 peak (Table S4, 
Figure S3). The deshielding effect of P-CD is commonly 
attributable to the formation of hydrogen bonding. These two 
results clearly reveal that dimethylamino and its adjacent hydroxyl 
(i.e., C3, C5, C6 and/or CI 2) of the molecule are involved in the 
inclusion by forming hydrogen bonding. 

The inclusion complex was visualized using molecular docking 
technique (Figure S4). At least four interactions (i.e., hydrogen 
bonding, hydrophobic interactions, van der Waals forces, and 
electrostatic interactions) are identified to form inclusion com- 
plexes but to different extents. Specifically, hydrogen bonding was 
observed to form between hydroxyl group/ ether bond at the wide 
rim of P-CD and C5/C12-hydroxyl groups of all the species with 
an exception of OTC ~. It supported the deshielding effect of CD 
in the NMR characterization. The number of hydrogen bonding 
involved foUowed an order of OTCH/ (2)>OTCH- 
(1) = OTC^" (l)>OTC" (0)) (Figure S4), which was in accordance 
with the decreasing order of Kq values (i.e., OTCH2^>"' 
OTCH->OTC^">OTC" (Table 1)). It can be concluded that 
hydrogen bonding plays the predominant role in the inclusion 
complexation of CD with the antibiotic. It was worthy noting that 
OTCH^ with less solubility and higher hydrophobicity [40] had 
stronger inclusion potential than OTC^~, though both species had 
the same number of hydrogen bonding. It was ascribed to the 
presence of hydrophobic interactions that commonly are recog- 
nized as driving forces for the inclusion [31]. Furthermore, total 
interaction energy of the complex, denoted a sum of van der Waals 
forces and electrostatic interactions, followed an order of 
OTC^">OTC">OTCH/>OTCH- that is independent of 
the R'c of the species (Table S5, Figure S5). This finding implied 
tiiat neither van der Waals forces nor electrostatic interactions 
played significant roles in the inclusion. 

Few attempts to investigate CD inclusion of lOCs with one- or 
two-pA;, also demonstrate that the molecular recognition of the 



species involves various intermolecular interactions depending on 
CD type and compound structure [41-44]. Sebestyen and cow- 
orkers summarize that the differences in the Kq of the species for 8 
amino acids are a consequence of the coordination of hydrogen 
bonding, steric effects and electrostatic interactions [42]. For 3 aro- 
matic amino acids, phenylalanine, tyrosine and tryptophan, the R'c, 
is in the range of 80-1 20 L/ mol for the anion while it is much smaller 
for the zwitterion and negligible for the protonated species. However, 
the differences were less pronounced for the other five amino acids. 

Adsorption of OTC onto CDPs 

The kinetics of adsorption was fast and the plateau reached within 
1 5-30 mill (Figure S6). The data of adsorption kinetics was well fitted 
to pseudo-second order model (Figure S7) relative to pseudo-first 
order model (Table S6). The rate controlling step was recognized as 
intraparticle diffusion for RMCDP and y-HP-CDP on the basis of the 
moving boundary model (Text SI) while liquid film diffusion for the 
other five polymers (Table S7). It is mainly due to rather high 
solubility [40] and thus hydrophiHc character of OTC, limiting mass 
transfer of the antibiotic to hydrophobic surfaces of materials. 

Figure SB represented adsorption isotherms of OTC. The 
adsorption constants derived from Langmuir [Ki) and FreuncUich 
{Ky) models were presented in Tables S8 and S9, respectively. In 
some cases, the sign of is negative (Table S8), indicating the 
poor application of Langmuir model. In contrast, FreuncUich 
model provided good fit to all data with i?^>0.92 (Table S9). The 
values of heterogeneity factor («) of FreuncUich model ranged from 
0.59 to 2.04 (Table S9). Obviously, most of them deviated from 
unit (n= 1) at which the FreuncUich isotherm translates to a linear 
relationship identical to the linear isotherm governed by hydro- 
phobic partitioning. Such deviation demonstrated that the OTC 
adsorption process was governed by complicated interactions 
other than hydrophobic interactions. Even for the same adsorbent, 
the values of n also varied with solution pH, indicating that 
adsorption interactions should differ at varying pH. The resulting 
Kp was of dimension, so that it would be a poor indicator for 
adsorption potential (Table S9). Instead of Aj, the dimensionless 
Aa, defined as the mean of at aU OTC levels (Text SI), was used 
to represent adsorption affinity of OTC in this study. Comparison 
of the values indicated that this term varied with pH and the 
maxima commonly approached at pH 5.5-7.5 (Figure 3). Here, 




H 1 . 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

3456789 10 
pH 

Figure 3. pH-dependent adsorption affinity of OTC onto CDPs. 

The scatter represents the experimental values at certain pH and the line 
represents the calculated ones at varying pH on the basis of Equation 2. 
doi:10.1371/journal.pone.0086228.g003 
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Table 2. Adsorption affinity {Kj, L/Kg) of the species onto CDPs. 







P-CDP 


RMCDP 


HPCDP 


y-CDP 


P-HP-CDP" 


P-Y-CDP" 


Y-HP-CDP" 


0TCH2* 


a 












1063 


OTCH- 


58 


246 


629 


111 


620 


141 


922 


OTC" 


27 


67 


178 


63 


137 


170 


342 


OTC^" ----- 



^approximately zero. 

"^The composite polymer was synthesized with an equimolar mixture of two CDs as double complexes. 
doi:l 0.1 371 /journal.pone.0086228.t002 



the values were comparable to tho.se reported for OTC onto 
other adsorbents, such as montmoriUonite, kaolinite and goethite 
(ca. 769, 100 and 900 L/Kg at pH 6.5-6.7) [45], and native 
montmoriUonite (ca. 1000 L/Kg at pH 5.0) [46]. 

Selective adsorption and immobilization of the species 

The term of the species (Table 2) was obtained through the 
LINGO optimization of the apparent Aj values of OTC at varying 
pH (Figure 3). There was good similarity between the apparent 
of OTC calculated on the basis of Equation 2 and obtained from 
batch adsorption experiments (Table SIO). OTCH^ commonly 
had higher values than OTC , with an exception of adsorbent 
P-y-CDP. By contrast, the values of OTCH2"*" and OTC"'" 
approached zero, except OTCH2^ onto y-HP-CDP with a 
value of 1063 L/Kg. 

FTIR characterization of the adsorbents (e.g., P-CDP) with 
OTC treatments at three pH values (i.e., pH 5.0, 7.0 and 9.0) 
further demonstrated selective adsorption and thus immobilization 
of the species onto the polymers (Figure S9). Specifically, 
stretching vibration band of hydroxyl of the polymer transferred 
from 3435 cm"' to 3385, 3415 and 3448 cm"' with OTC 
treatments at pH 5.0, 7.0 and 9.0, respectively. Stretching band of 
C-OH at 1034 cm"' was widened for OTC treatments at pH 5.0 
and 7.0, while it disappeared for that at pH 9.0. Moreover, C-O-C 
stretching vibration shifted from 1154 cm" to 1158 cm" 
(pH 5.0), 1160 cm"' (pH 7.0) and 1165 cm"' (pH 9.0), respec- 
tively. The species-specific adsorption may be attributed to both 
ionization of OTC at varying pH and different adsorption 
interactions of the species, similar to the finding of the pH- 
dependent interactions of OTC in clay and organic matter [46] . 

In a recent study [47], surface coverage of the species has been 
highlighted for citric acid onto Ti02 anatase nanoparticles, and 
the fully deprotonated species is preferred to immobilization. In 
this study, immobilization potential (i.e., adsorption amounts) of 
the species for OTC onto the polymers was calculated on the basis 
of the and aqueous fractions of the species. The results 
demonstrated the immobilization of the species was a function of 
pH (Figure 4). For adsorbent y-HP-CDP, the species immobilized 
was OTCHz"" (9.7-68.2%) and OTCH" (31.8-90.2%) at pH 3.0- 
4.3, OTCH- (91.8-98.8%) atpH 4.4-6.7, OTCH" (18.2-89.9%) 
and OTC" (10.0-81.8%) at pH 6.8-8.5, and OTC" (84.4- 
95.3%) at pH 8.6-10.0. For the other six adsorbents onto which 
immobilization of OTCHj^ was excluded, the species immobilized 
was OTCH- (90.3-100%) at pH 3.0-6.2, OTCH" (2.2-97.9%) 
and OTC" (2.1-97.8%) at pH 6.3-9.3, and OTC" (90.0-98.5%) 
at pH 9.4-10.0. Obviously, the difiFerences in immobilization of 
the species are responsible for the apparent pH-dependent 
adsorption of OTC. 



Recognition of interactions of the species of OTC onto 
CDPs 

As it is well recognized, CD inclusion is the dominant 
adsorption interaction for CD-based polymers toward most guest 
compounds [17-19,38]. In most cases, high CD content corre- 
sponds to high adsorption affinity [17-19,38]. In contradiction to 
this finding, the of the species in this study (Table 2) was not 
correlated with CD content of the polymers that is compiled in 
Table 1 in our previous work [19]. Especially, y-HP-CDP that has 
the second lowest CD content [19] exhibited the highest affinity to 
the species OTCH^^, OTCH*, and OTC" (Table 2). Such 
deviation suggested that the adsorption of the species be governed 
by a combination of CD inclusion and other mechanisms. 

The adsorption interactions were attempted to distinguish using 
multiple linear regression (MLR) analysis for the of OTCH^ 
and OTCH" with CDPs properties. CD content, cross-linking 
degree, swelling ratio, particle size, pore size, surface area, pore 
volume, and species-specific inclusion constant were respectively 
set as variables xi, X2, X3, X4, X5, xg, x-j, and xg. The MLR analysis 
was performed by the stepwise method in the SPSS Statistic 
program (SPSS 17.0 version) (Equations 3 and 4). CoUinearity 
between variables was excluded in terms of variance inflation 
factor <10. The analysis was repeated (Equations 4 and 5) when 
substituting variables xi and Xjj with variable xg (denoted xi xx^). 
Variable Xg, relating to both CD content and species-specific 
inclusion constant, may be indicative of effective inclusion. 




-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 

3456789 10 



Figure 4. Species-specific fractions of OTC immobilized onto 
CDPs. The hollow with a cross, solid, and hollow scatter represent the 
fractions of OTCH2*, OTCH* and OTC", respectively. 
doi:10.1371/journal.pone.0086228.g004 
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OTCH± : i:d=0.879.x:7 (i?^ =0.765) (3) 
OTC : iiCd =0.819x2 (i?^ =0.670) (4) 

OTCH± : ^d = l-075.T7 + 0.453.x:9 {R^ = OmO) (5) 

The adsorption of OTCH^ was correlated with pore volume 
and effective inclusion capacity (Equations 3 and 5). It highlights 
the presence of pore-filling [48,49] and CD inclusion mechanisms 
for the species OTCH^. Comparison of standardized coefficients 
showed that pore-fUling mechanism contributed more than CD 
inclusion one. It should be noted that the pore of the polymer is 
composed of pore channel and CD cavity. Analysis of the complex 
of the species with CD (Figures S4, S5) or EPI oligomers (Figure 



SIO) demonstrated that both units could bind the species via 
hydrogen bonding and/or a-Tt interactions. The substitution of 
variable Xg rendered a significant increase in the regression 
coefficient {R ) (Equation 5 versus Equation 3). This improvement 
suggested that variable xg was a good indicator of inclusion 
potential relative to variables xi and xg. A concomitant of the 
regression improvement was an increase in the weight coefficient 
of pore volume. It is probably due to the dual roles (i.e., inclusion 
complexation and pore donator) of CD, both which are derived 
from the cavity of CD. 

The adsorption of OTC was mainly related to cross-linking 
degree (Equation 4). High cross-linking degree commonly causes a 
dense network structure [19,50]. Hydrogen bonding forces and a- 
71 interactions were identified for the species OTC binding to EPI 
oligomers (Figure SIO). This correlation probably illustrates the 
presence of network capture mechanism. In contrast, the 
substitution of variable xg didn't affect the MLR analysis of 
OTC . This indicated an insignificant role of CD inclusion for 
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Figure 5. Multiple linear regression analysis of apparent adsorption of OTC at varying pH with CDPs properties. I, both CD content 
and inclusion constant were set as independent variables; and II, CD contentxinclusion constant was set as an independent variable. 
doi:1 0.1 371/journal.pone.0086228.g005 
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Hydrogen bonding and van der Waals interactions 


Hydrophobic interactions 
Electrostatic attractions , , ct-ji forces , 


^"^^^^pH 3-4 
' ^^^^^^^^^^^^ 




OTCH2+ OTCH± 


OTC OTC^ 



Figure 6. Species-specific interactions involved in the adsorption of OTC at varying pH onto CDPs. 

doi:10.1371/journal.pone.0086228.g006 



OTC , which was in agreement with the finding of its poor CD 
inclusion potential (Table 1). 

For the other two species that are present in acidic or basic 
solutions (i.e., OTCH2^ and OTC^ ), the MLR analysis wasn't 
performed because the almost approached zero (Table 2). The 
poor adsorption affinity of the two species could be attributed to 
the fact that the adsorption sites (i.e., CD cavity) were not readily 
accessible for the two species, because they were difficult to diffuse 
to hydrophobic surfaces (Table S7). Furthermore, they could 
interact with hydroxyl groups in the network but in the less 
magnitude, as compared with the other two species. It is due to the 
decrease in the number of hydrogen bonding for the species 
following the order of OTCH" (10)>OTCH- (8)>both species 
(3) (Figure SIO). 

The nature of pH-dependent adsorption of OTC 

The MLR analysis of apparent of OTC with CDP properties 
was performed, in which the apparent values of OTC was 
calculated on the basis of Equation 2 (showed as lines in Figure 3) 
when pH varied from 3.0 to 10.0 with an interval of 0.2. The 
apparent adsorption of OTC at varying pH was correlated with 
one to five of all eight factors of concern, when CD content and 
inclusion constant were set as independent variable (Figure 51). 
Surprisingly, the substitution of variable reduced the number of 
influencing factors, only including pore volume and/or effective 
inclusion, while regression coefficients {R^) were comparable 
(Figure 511). 

A combination of the MLR analysis for OTC (Figure 511) and 
molecular recognition in adsorption of the species (Equations 3-5) 
may illustrate pH-dependent adsorption of OTC onto CDPs 
(Figure 6). Specifically, the term correlates with pore volume at 
pH 3.0-4.2, where OTCH* and OTCHj"^ prevail and are 
adsorbed via pore-filling mechanism. As the solution pH rises up 
to 6.8, pore volume and effective inclusion capacity become the 
influencing factors corresponding to pore-fiUing and CD inclusion 
mechanisms, respectively. It is reasonable because OTCH^ is the 
predominant species in the range of pH. However, the standard- 
ized coefficients of both influencing factors (Figure 511) gradually 
deviate from those in Equation 5 with pH increasing. It is probably 
due to an increase in the fraction of OTC over pH. Surprisingly, 
the term pore volume is the sole influencing factor at pH 7.0-10.0, 



indicating the presence of pore-fiUing mechanism. It is in 
contradiction to the fact that the species OTC~ predominates at 
pH 8.0-9.0 and its adsorption proceeds via network capture 
mechanism. This may be due to a combination of pore-flUing, CD 
inclusion and network capture mechanisms, associated with the 
coexistence of two or three species in this pH range. 

Selection of adsorbent to separate OTC from water 

The selection of the polymers may be determined on the basis of 
the understanding of the nature of pH-dependent adsorption of 
OTC (Figure 511). Adsorbent y-HP-CDP with advantages of high 
pore volume and efficient inclusion potential was selected to 
separate and remove OTC from simulated wastewater (pH 5.0) 
and natural water (pH 7.0) (Text S4). The adsorbent removed 
80% of OTC at 20 mg/L (pH 5.0) (Figure 7). In contrast, other 
conventional natural and synthetic adsorbents are observed to 
exhibit similar or lower adsorption efficiency for OTC, such as 
activated carbon (ca. 68%) [7,51], multi-wall carbon nanotubes 
(ca. 70%) [52], and aluminum oxide (44%) [53]. Furthermore, the 



pH=7.0 lOng/L 
pH=5.0 20mg/L 




1 2 3 4 5 

Repetition times 



Figure 7. Removal efficiency of OTC by adsorbent y-HP-CDP. 

doi:10.1371/journal.pone.0086228.g007 
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Table 3. Effect of coexisting components on removal 
efficiency of OTC by y-HP-CDP at pH 5.0. 




Removal efficiency, % 


OTC, 20 mg/L (control) 


80.7 ±0.9 


Ca^+, 1 50 mg/L 


80.5 ±0.2 


Ca^+, 50 mg/L 


81.3±0.4 


Mg^+ 50 mg/L 


80.4 ±0.1 


Mg^+, 8 mg/L 


80.8±0.8 


DOM, 10 mg/L 


80.7 ±1.0 



doi:l 0.1 371 /journal.pone.0086228.t003 



polymer exhibits a removal efficiency of above 70% for the 
antibiotic even at environmentally relevant levels (Figure 7). 

In addition, ubiquitously existing Ca^"^ (e.g., at 150 and 50 mg/ 
L), Mg^"^ (e.g., at 50 and 8 mg/L) and DOM (e.g., at 10 mg/L) 
that are natural water constituents and usually are indicative of 
environmental matrix effects did not deteriorate removal efficiency 
of OTC in wastewater (pH = 5.0) (Table 3). This is consistent with 
the findings of other studies about the removal of neutral 
pollutants by cyclodextrin materials [19,54]. In contrast, many 
works have showed negative effects of cations (Na"^, Ca^"^, Mg^"*") 
and DOM on removal rate of tetracyclines by conventional 
adsorbents (e.g., Na-montmoriUonite [55], humic acids [56], 
graphene oxide [57], CNTs [58], kaolinite [59], and clays [46]). 
In some cases, the decrease in removal rate approaches about 50- 
90% [58,59]. 

Conclusions 

The finding of the species-specific adsorption reveals the pH- 
dependent adsorption behavior of antibiotic oxytetracycline onto 
adsorbent CDPs. The adsorption sites of the polymer such as CD 
cavity, pore channel, and network are accessible for the species to 
different extents, associated with the coordination of a variety of 
weak interactions (e.g., hydrogen bonding, hydrophobic interac- 
tions, van der Waals forces and electrostatic interactions). The 
species OTCH^ commonly has high adsorption affinity via pore- 
fiUing and CD inclusion mechanisms, OTC exhibits moderate 
affinity via network capture mechanism and the other two species 
are poorly adsorbed due to the limiting accessibility of the 
inclusion sites as well as other adsorption sites. The resulting 
immobilization of the species onto CDPs is determined by the 
and fraction of the species in solution, representing a function of 
pH. This may be the cause for the pH-dependence of adsorption 
of OTC. Moreover, the mathematical equations are derived to 
quantitatively relate adsorption affinity of OTC at varying pH to 
adsorbent properties. These relationships, together with the mole- 
cular recognition into adsorption of the species, provide theoretical 
guidelines for the selection of CDPs to remove the antibiotic. 
Additionally, this study indicates that CDPs are promising 
adsorbents for lOCs removal, taking advantage of discernible 
adsorption interactions, high removal efficiency, and strong 
resistance to the interference of coexisting natural components. 

Supporting information 

Text SI Adsorption kinetics and isotherms of OTC. 

(DOC) 

Text S2 Determination of the stability of OTC. 

(DOC) 



Text S3 LINGO optimization program of species-spe- 
cific interaction constant. 

(DOC) 

Text S4 Removal of OTC from simulated waters. 

(DOC) 

Table SI Related information of cyclodextrin. 

(DOC) 

Table S2 Inclusion complexation of CD with OTC at 
varying pH. 

(DOC) 

Table S3 Linear correlation between Kc values calcu- 
lated and experimental. 

(DOC) 

Table S4 chemical shifts (ppm) of OTC, CD and CD- 
OTC protons. 

(DOC) 

Table S5 Intermolecular interactions of CD and the 
species. 

(DOC) 

Table S6 Parameters of the pseudo -first and -second 
order model. 

(DOC) 

Table S7 Correlation coefficients (R^) of the moving 
boundary model. 

(DOC) 

Table S8 Adsorption parameters of OTC fitted to 
Langmuir model. 

(DOC) 

Table S9 Adsorption parameters of OTC fitted to 
Freunillich model. 

(DOC) 

Table SIO Linear correlation between values calcu- 
lated and experimental. 

(DOC) 

Figure SI UV-Vis spectra of CD-OTC complexes at 
varying pH. 

(TIF) 

Figure S2 FTIR spectra of CD, OTC and CD-OTC 
complexes. 

(TIF) 

Figure S3 *H NMR spectra of OTC, CD and CD-OTC 
complexes (dissolved in D2O). HOD at 4.617 ppm was taken 
as a reference. 
(TIF) 

Figure S4 Molecular docking of CD with the species. 

The green dotted line represents hydrogen bonding. 
(TIF) 

Figure S5 Qualitative comparison of intermolecular 
interactions between CD and the species. The blue, green, 
and red isosurfaces are indicative of hydrogen bonding, van der 
Waals forces, and repulsion interactions, respectively. 
(TIF) 

Figure S6 Adsorption kinetics of OTC onto CDPs at 
pH 7.0. 

(TIF) 
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Figure S7 Pseudo-second order model for adsorption of 
OTC onto CDPs (pH 7.0). 

(TIF) 

Figure S8 Adsorption isotherms of OTC onto CDPs at 
varying pH. 

(TIF) 

Figure S9 FTIR spectra of CDP adsorbing OTC at 
respective pH 5.0, 7.0 and 9.0. 

(TIF) 

Figure SIO Molecular docking of the EPI oligomer with 
the species. The green dotted line represents hydrogen bonding 



and the orange line represents cy-n interaction. The numbers in 
parentheses represent ring-opened, tricyclo-, pentacyclo-, and 
heptacyclo-forms, respectively. 
(TIF) 
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